Equilibrium Distribution of the Inherent States and 
their Dynamics in Glassy Systems and Granular Media 
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The present paper proposes a Statistical Mechanics ap- 
proach to the inherent states of glassy systems and granular 
materials, following the original ideas developed by Edwards 
for granular materials. Two lattice models, a diluted Spin 
Glass and a system of hard-spheres under gravity, introduced 
in the context of glassy systems and granular materials, are 
evolved using a "tap dynamics" analogous to that of exper- 
iments on granular materials. The asymptotic macrostates, 
reached by the system, are shown to be described by a single 
thermodynamical parameter, and this parameter to coincide 
with the temperature, called the "configurational tempera- 
ture", predicted assuming that the distribution among the 
inherent states satisfies the principle of maximum entropy. 

The thermodynamics of macroscopic systems evolving 
at equilibrium is well described by Statistical Mechan- 
ics. However there are many systems, typically found in 
"frozen states", where they do not evolve at all. These 
are, for example, supercooled liquids quenched at zero 
temperature in states, called inherent states [1,2], corre- 
sponding to the local minima of the potential energy in 
the 3N-dimensional configuration space of particle coor- 
dinates. Granular materials [3] at rest are another im- 
portant example of system frozen in mechanically stable 
microstates. Grains are "frozen" because, due to their 
large masses [3], the thermal kinetic energy is negligible 
compared to the gravitational energy; thus the external 
bath temperature, T{,ath, can be considered equal to zero 
(by analogy with supercooled liquids, we call these me- 
chanically stable configurations inherent states). 

In this paper, following the original ideas by Edwards 
for granular materials [4] we attempt to develop a unified 
Statistical Mechanics approach for the inherent states of 
glassy systems and granular materials along the line of 
Ref. [5] . The connection between Edwards approach and 
recent developments on glass theory has received much 
attention [6-12]. 

The first step is to introduce a suitable dynamics which 
allows to explore the configurations of the inherent states. 
In granular materials the dynamics, from one stable mi- 
crostate to another, can be induced by sequences of 
"taps" , in which the energy is pumped into the system 
in pulses. Due to inelastic collisions the kinetic energy is 
totally dissipated after each tap, and the system is again 
frozen in one of its inherent states [13]. Similarly, in glass 
formers at zero temperature the dynamics, from one in- 
herent state to another, can be induced by sequences of 
taps, where each tap consists in raising the bath temper- 
ature and, after a lapse of time tq, quenching it back to 



zero. By repeating the process cyclically the system ex- 
plores the space of the inherent states [5,7-9,14,15]. For 
a tap of infinite length (tq —> cxd) the way to explore 
the inherent states coincides with the one used in [1,2] 
for a system of Lennard Jones mixture. In the approach 
of Barrat et al. [6] the system instead evolves in an out 
of equilibrium quasi-stationary state at an external very 
low bath temperature. In the limit of zero external tem- 
perature the system explores the inherent states. 

Here we consider a diluted Spin Glass and a system of 
hard-spheres under gravity, introduced in the context of 
glassy systems and granular materials, which are evolved 
using a "tap dynamics" . We show that the systems reach 
a stationary or quasi-stationary state which can be char- 
acterized by a single thermodynamical parameter, de- 
fined through the static fluctuation-dissipation relation, 
and this parameter coincides with the "configurational 
temperature" , predicted by the Edwards' hypothesis of 
a flat measure for the microstate distribution. We also 
show that time averages over the dynamics can be re- 
placed by ensemble averages over such measure. 

We first consider the Frustrated Lattice Gas model in 
three dimensions {3D). The model was recently intro- 
duced to describe glass formers [16] and, in presence of 
gravity, granular materials [15,17,18]. It is made of parti- 
cles with a twofold orientation (i.e., pointing in only two 
allowed directions), displaced on a cubic lattice (of linear 
size L = 8 and overall density p = J^i ^i/L^ — 0.65 and 
0.75), and interacting via a quenched potential: 



H = J ^ fij {Si, Sj)ninj , 



(1) 



where rii — 0,1 whether site i is empty or filled by a 
particle. Si — ±1 is a variable associated to the particle 
orientation, J the amplitude of the interaction poten- 
tial between neighbours. The shape factor fij{Si,Sj) = 
1 — CijSiSj (where e.y = ±1 are quenched and random 
variables) is or 1 depending whether the relative ori- 
entations {Si, Sj) are favoured or not when neighbouring 
sites i and j are both occupied (for a review see [17]). 

We also consider a 3D system of hard-spheres subject 
to gravity, where the centers of mass of grains are con- 
strained to move on the sites of a cubic lattice, as de- 
picted in the upper inset of Fig. 4. Its Hamiltonian is 
given by eq.(l), plus a gravitational term, with /^ = 1 
if i and j are nearest neighbours, fij = elsewhere and 
J — > cxD (in this case grains have no orientation and Si is 
redundant \/i). 



In both models the value of particle density is fixed, 
and a Monte Carlo tap dynamics, which allows the sys- 
tem to explore its inherent states, is applied. During the 
dynamics, the system cyclically evolves for a time tq (the 
tap duration [19]) at a finite value of the bath temper- 
ature, Tr (the tap amplitude), and it is suddenly frozen 
at zero temperature in one of its inherent states (at zero 
temperature the system does not evolve anymore if the 
energy cannot be decreased by one single particle move- 
ment). After each tap, when the system is at rest, we 
record the quantities of interest. The time, t, considered 
is therefore discrete and coincides with the number of 
taps. 

Let's first discuss the results about the Frustrated Lat- 
tice Gas model for density, p — 0.65, since very similar 
features are found for p = 0.75 and in the hard-sphere 
system (described later on). 

Interestingly, under the tap dynamics the system 
reaches a stationary state (see Fig.l) for all the values 
of Tr (and tq) we considered. During the tap dynamics, 
in the stationary state, we have calculated the time av- 
erage of the energy, E, and its fluctuations, AE"^. We 
show the results in Fig. 2 as function the tap amplitude, 
Tr (for several values of the tap duration, tq). Appar- 
ently, Tr is not the right thermodynamical parameter, 
since sequences of "taps" with different tq give differ- 
ent values of E{Tr, tq) and AE'^{Tr, tq). However, if the 
stationary states corresponding to different tap dynamics 
(i.e., different Tr and tq), are indeed characterised by a 
single thermodynamical parameter, all the curves should 
collapse onto a single master function when AE^ is para- 
mctrically plotted as function of E. This data collapse 
is in fact found and shown in the lower inset of Fig. 3. 
This is a prediction which could be easily checked in real 
granular materials (where one could consider the density 
which is easier to measure than the energy). 

In the Frustrated Lattice Gas for density, p = 0.75, as 
much as in the hard-sphere model, for low values of the 
tap amplitude, Tr, the system reaches a quasi-stationary 
state, where one time quantities decay as the logarithm of 
time. In this case the average over the time is performed 
over a time interval of the tap dynamics such that the 
energy is practically constant. By performing then the 
same procedure described in the stationary case, we find 
(see Ref. [20]) again a collapse of data as in Fig. 3. 

The thermodynamical parameter, Pfd, is defined 
apart from an integration constant, /3o, through the 
fluctuation-dissipation relation: 



dE 



AE^. 



(2) 



By integrating eq. (2), /3/d — /?o can be expressed as 
function of E, where the integration constant Po can be 
determined independently [20]. In Fig. 3, E as function 
of Tfd = PJd is shown. The corresponding results for 
the hard-sphere system are shown in Fig. 4. In the last 



model we have also checked that the system density on 
the bottom layer, pt, and the density self-overlap func- 
tion, Q depend only on f3fd (see Ref. [20]), confirming 
that a unique thermodynamical parameter is enough to 
describe the system macrostates. 

We have found that the fluctuations of the energy in 
the stationary state depend only on the energy, E, and 
not on the past history. If all macroscopic quantities 
depend only on the energy, E, or on its conjugate ther- 
modynamical parameter, (3fd, the stationary state can be 
genuinely considered a "thermodynamical state" . If this 
is the case one can attempt to construct an equilibrium 
statistical mechanics, as originally suggested by Edwards 

[4]. 

More precisely we ask in the stationary regime what 
is the probability distribution, Pr, of finding the system 
in the inherent state r of energy Er (see [5]). We as- 
sume that the distribution is given by the principle of 
maximum entropy, S = —^^ PrlnPj., under the condi- 
tion that the average energy is fixed: E = ^^PrEr- 
Thus, we have to maximise the following functional: 

I[Pr] - - E. PrlnPr " (^confiE - E. ^r^r)- Here Pconf 

is a Lagrange multiplier determined by the constraint on 
the energy and takes the name of "inverse configurational 
temperature". This procedure leads to the Gibbs result: 



Tr- 



-/3c„„/S. 



(3) 



where Z — Er e~'^'=°"^^''. Using standard Statistical Me- 
chanics it is easy to show that, in the thermodynamic 
limit, the entropy S and Pconf are also given by: 



S^lnQ{E), Pconf 



dlnft 
dE 



(4) 



where fl{E) is the number of inherent states correspond- 
ing to energy E. 

It is possible to show [5] that in the particular case in 
which the particles density p is constant, Tconf is simply 

related to the "compactivity" , X — (w) , introduced 
by Edwards in his seminal papers [4]. In general X and 
Tconf are two independent variables. 

If the distribution in the stationary state coincides 
with eq.(3) the time average of the energy, E{Pfd), 
recorded during the tap dynamics, must coincide with 
the ensemble average, {E){f3conf), over the distribution 
eq.(3). To calculate the average (E), as function of 
(3conf, we have introduced an auxiliary hamiltonian (see 
also [6]) H'{{S„n^}) = H{{S,,n,}) + S{{S^,n,}), where 
7i({5'i,ni}) is the hamiltonian (1), and (5({S'i,ni}), is 
zero, if the configuration is an inherent one, and infinite, 
otherwise. In this way the canonical distribution for this 
Hamiltonian gives a weight, e~^''°"f''^ , which is equal to 
g-Pcon/H^ for the inherent configurations, and zero oth- 
erwise, reproducing the distribution cq.(3). With this 
auxiliary hamiltonian, using standard Monte Carlo sim- 
ulations, we have then calculated {E){(3conf)- Fig.s 3, 4 



outline a very good agreement between {E){Pconf) and 
E{(ifd) (notice that there are no adjustable parameters). 
The same agreement is found in the Frustrated Lattice 
Gas model for p = 0.75 and in the hard-sphere system for 
the other quoted observables, pb, and Q (see Ref. [20]). 

In the approach of Ref. [6] the system explores the 
inherent states evolving in an out of equilibrium quasi- 
stationary state at a very low bath temperature; in this 
approach the configurational temperature is expected to 
coincide with the "dynamical temperature", Tdyn, which 
appears in the extension of the fluctuation-dissipation re- 
lation to the out-equilibrium case. One of the differences 
with the approach used here is that using the tap dynam- 
ics it is also possible to explore low density inherent states 
in a stationary regime and not only the off-equilibrium 
"glassy regime" . For istance, the Frustrated Lattice Gas 
model at density p ~ 0.65 (one of the cases here studied) 
is never found in an out of equilibrium quasi-stationary 
state (at any finite value of the bath temperature the 
system quickly reaches the equilibrium state). 

In conclusion, in the context of models for glasses and 
granular materials, we have obtained two different re- 
sults. First, we have shown that the stationary states 
reached by the system under the tap dynamics among the 
inherent states are not dependent on the past history and 
can be considered as a thermodynamical state character- 
ized by a single thermodynamical parameter, Tfd, defined 
through the fluctuation-dissipation relation. Second, Tfd 
coincides with Tconf, predicted assuming that the distri- 
bution among the inherent states satisfies the principle 
of maximum entropy under the constraint that energy 
is fixed. Moreover ensemble average coincides with time 
average over the tap dynamics. In particular we have 
found that, by using Tconf as a state parameter, the ob- 
servables recorded in different tap sequences (different 
amplitude and duration of taps) fall onto universal mas- 
ter curves, and the curves turn out to coincide with the 
ones predicted by the distribution eq.(3). 



[4] S.F. Edwards and R.B.S. Oakeshott, Physica A 157, 
1080 (1989). A. Mehta and S.F. Edwards, Physica A 157, 
1091 (1989). S.F. Edwards, in "Disorder in Condensed 
Matter Physics" page. 148, Oxford Science Pubs (1991); 
and in Granular Matter: an interdisciplinary approach, 
(Springer- Verlag, New York, 1994), A. Mehta ed. 
[5] A. Coniglio and M. Nicodemi Physica A 296, 451 (2001). 
[6] A. Barrat, J. Kurchan, V. Loreto and M. Sellitto, Phys. 
Rev. Lett. 85, 5034 (2000). Phys. Rev. E 63, 051301 
(2001). 
[7] J.J. Brey, A. Prados, B. Sanchez-Rey, Physica A 275, 
310 (2000); Physica A 284, 277 (2000). A. Prados and 
J.J. Brey, cond-mat/0106236. 
[8] A. Lefevre, D. S. Dean, Phys. Rev. Lett. 86, 5639 (2001); 

cond-mat/0106220. 
[9] J. Berg, S. Franz and M. Sellitto, cond-mat/0111485. 

[10] L. Berthier, L.F. Cugliandolo and J.L. Iguain, Phys. Rev. 
E 63, 051302 (2001). 

[11] H. A. Makse and J. Kurchan, Nature 415, 614 (2002). 

[12] V. Colizza, A. Barrat and V. Loreto cond-mat/0 111458. 

[13] J.B. Knight et al., Phys. Rev. E 51, 3957 (1995). 

[14] J. Berg, A. Mehta, cond-mat/0012416. 

[15] M. Nicodemi, A. Coniglio, H.J. Herrmann, Phys. Rev. E 
55, 3962 (1997); J. Phys. A 30, L379 (1997). A. Coniglio 
and H.J. Herrmann, Physica A 225, 1 (1996). 

[16] A. Coniglio, A. de Candia, A. Fierro, M. Nicodemi, Jour. 
Phys.; Cond. Mat. 11, A167 (1999). 

[17] A. Coniglio, M. Nicodemi, J. Phys.: Cond. Matt. 12, 
6601 (2000); Phys. Rev. Lett. 82, 916 (1999). M. 
Nicodemi, Phys. Rev. Lett. 82, 3734 (1999). 

[18] J.J. Arenzon, J. Phys. A. 32, L107 (1999). 

[19] To is measured in Monte Carlo steps (MCS), where 
1 MCS corresponds to A'^ attempts to move a particle 
randomly chosen, and A'' is the number of particles. 

[20] A. Fierro, M. Nicodemi and A. Coniglio, in preparation. 



We thank S.F. Edwards for useful discussions. 
This work was partially supported by the TMR- 
ERBFMRXCT980183, INFM-PRA(HOP), MURST- 
PRIN 2000. The allocation of computer resources from 
INFM Progetto Calcolo Parallelo is acknowledged. 



[1] F.H. Stillinger T.A. Weber, Phys. Rev. A 25, 978 (1982). 

S. Sastry, P.G. Debenedetti, F.H. Stillinger, Nature 393, 

554 (1998). 
[2] F. Sciortino, W. Kob, P. Tartagha, Phys. Rev. Lett. 83, 

3214 (1999). F. Sciortino and P. Tartagha, Phys. Rev. 

Lett. 86, 107 (2001). 
[3] H.M. Jaeger, S.R. Nagel and R.P. Behringer, Rev. Mod. 

Phys. 68, 1259 (1996). 



1.0 
0.8 
3 0.6 
[?"0.4 
0.2 
0.0 



10^ 
10' 


"■■^ 










80 

^ ©60 

\ 20 


\ 




• 
• 

\ 

V 




3° 10 


10 


'lo'io^io 
t 








n t„-io- 

'^ ' t -10"' 


10° 10' 
Tr 


• t.-io' 




V-i-i, 


H^ 


•...-. 



10° lO' 10^ 10^ lO" 10^ 

t-t„ 



FIG. 1. Upper inset The energy, E{t), of inherent states 
as a function of the "tap" number t, during a tap dy- 
namics with a tap amplitude Tr = 1.25 J and tap dura- 
tion To = 1 MCS, in the Frustrated Lattice Gas of the 
text. Main frame The self-scattering two times function, 



F^{t,t^) = J2i(^' 



<i-lPi(-t)-fi(t^)] 



I pL as a function of the num- 



ber t - tu, (7r = 1 J, TO = 1 MCS and q = 7r/4). Lower 
Inset The equilibration time, r, as a function of Tr (for 
To = 1 MCS): T diverges at low Tr. After a transient, 
E{t) reaches its time independent asymptotic value, and 
Fq = Fq{t — tw) depends only on the difference of i — t^, 
showing that the system has reached a stationary state (our 
data are averaged up to 32 noise realizations). 
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FIG. 2. The time average of the energy, E, and (inset) 
its fluctuations, AE'^ , recorded at stationarity during a tap 
dynamics with tap amplitude, Tr, in the Frustrated Lattice 
Gas model. The four different curves correspond to sequences 
of tap with different values of the duration of each single tap, 
TO = 100, 10, 5, 1 MCS (from top to bottom). 
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FIG. 4. Main frame The time average, E, and the ensem- 
ble average over the distribution eq.(3), {E) (black filled dia- 
monds), plotted respectively as a function of Tjd and Tconf, 
in the 3D hard-sphere system under gravity described in the 
text (and schematically depicted in the upper inset). Lower 
Inset The temperature Tfd = PJ^ defined by eq.(2) as func- 
tion of Tr for to = 500, 10, 5 MCS (from top to bottom). The 
straight line is the function Tfd = Tr- 
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FIG. 3. Lower Inset The time averages of energy fiuctu- 
ations, AE^, when plotted as function of the time average 
of energy, E, collapse on a single master function for all the 
different values of tap amplitude and duration, Tr and to, 
shown in Fig. 2. Main frame The time average, E, and the 
ensemble average over the distribution eq.(3), (E) (black filled 
diamonds), plotted respectively as a function of Tfd = PJ^ 
and Tconf, in the Frustrated Lattice Gas model. Upper In- 
set The temperature Tfd = [37^ defined by eq.(2) as function 
of Tr for to = 100, 10, 5, 1 MCS (from top to bottom). The 
straight line is the function Tfd = Tr. 



